clc;clear all;close all;

R = 120;
C = 33e-12;
tau = R*C;

a = [tau 1];
b = [1];

Fmax = 100e6;
Fs = 1000;
f = 0:Fs:Fmax;

%ачх и фчх

H = freqs(b,a,2*pi*f);

figure(1);subplot(2,1,1);
plot(f/1e6,20*log10(abs(H)));
grid on;
ylabel('|H|');
xlabel('f');
subplot(2,1,2);
plot(f/1e6,rad2deg(unwrap(angle(H))));
ylabel('\phi');
xlabel('f');
pause(1);
%импульсная характеристика
sys = tf(b,a);
[y,t] = impulse(sys);
figure(2); plot(t, y);
xlabel('t, s'); ylabel('h(t)');
grid on

%нули и полюсы
[z,p,k]=tf2zp(b,a);
figure(3);
zplane(z,p);

%дискретная импульсная характеристика
[bd,ad] = bilinear(b,a,8*10e6);

h_imp = impz(bd,ad,15);
figure(4);stem(h_imp);
xlabel('k');ylabel('h_imp');grid on;

%свертка
x = [0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0];
y = filter(bd,ad,x);
figure(5)
stem(1:length(y), y); hold on
stem(1:length(x), x, 'r'); hold off
grid on; legend('y', 'x'); xlabel('k');


%нули и полюсы
[z,p,k]=tf2zp(bd,ad);
figure(6);
zplane(z,p);

%фильтр
figure(7);
t = 0:1/10:100*Fs;
x = sin(2*pi*t*40e6)
y = filter(bd,ad,x);
stem(1:length(y), y); hold on
stem(1:length(x), x, 'r'); hold off
grid on; legend('y', 'x'); xlabel('k');